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We study a multiple-spin exchange model on a triangular lattice, which is a possible model for 
low-density solid ^He films. Due to strong competitions between ferromagnetic three-spin exchange 
and antiferromagnetic four-spin one, the ground states are highly degenerate in the classical limit. 
At least 2^''^-fold degeneracy exists on the L x L triangular lattice except for the SO(3) symmetry. 
In the magnetization process, we found a plateau at m/msat ~ 1/2, in which the ground state is 
uuud state (a coUinear state with four sublattices). The 1/2-plateau appears due to the strong four- 
spin exchange interaction. This plateau survives against both quantum and thermal fluctuations. 
Under a magnetic field which realizes the uuud ordered state, a phase transition occurs at a finite 
temperature due to breakdown of translational symmetry. We predict that low-density solid ^He 
thin films may show the 1/2-plateau in the magnetization process. Experimental observation of the 
plateau will verify strength of the four-spin exchange. It is also discussed that this magnetization 
plateau can be understood as an insulating-conducting transition in particle picture. 



I. INTRODUCTION 

In localized fermion systems, the magnetic interac- 
tion comes from permutations of particles.El'B Multiple- 
spin exchanges, e.g. cyclic exchanges of three or four 
spins, have been revealed to be strong in nuclear mag- 
netism of tW|Q-dimensional (2D) i^olid "^He films. Many 
experimentalBia and thcorcticalErb2l studies suggested 
that exchange interactions of more-than-two spins are 
dominant in this system. A ferromagnetic behavior 
changes to antLEecromagnetic one when the coverage of 
^He decreases .E1It13 This tendency can be understood in 
terms of multiple-spin exchange (MSE): in fuUy packed 
systems, the three-spin exchange is depiinantQ and it is 
ferromagnetic (as shown by Thoulesscl), and, in loosely 
packed systems, the four- and six-spin exchanges be- 
come strong and favor antiferromagnetism. Effects of 
multiple-spin exchanges are not yet-iully understood es- 
pecially for the low-density region. E£IeI For example, re- 
cent specific-heat data at low densities show a peculiar 
behavior, which have a double-peak structure, and they 
also show that the ground state seems to be spin liquid 
(disordered) ff^d the spin excitation gap is vanishing (or 
quite small) .eI 

The multiple-spin exchange (MSE) model has been 
studied to describe thcj-iiuclear magnetism of three- 
dimensional solid ■^HeJlH A general fosfft of the 
spin Hamiltonian of quantum solid isElfl Ti, = 
~I]n(~l)"'^"I]p„ -P"' where P„ and J„(< 0) denote 
cyclic permutation of n spins and its exchange con- 
stant, respeptpely. For the 2D system, recent theoretical 
calculationsBiij and experimental measurementsu found 
that the exchange frequencies satisfy | J3 1 > | J2 1 > | J4 1 ^ 
I "^6 1 ^ I "^5 1 on the triangular lattice. In this paper, we 
consider a spin model with the two-, three-, and four-spin 
exchanges on the triangular lattice, which is the simplest 
2D MSE model. Since the three-spin exchange can be 



transformed to the two-spin ones, the Hamiltonian can 
be written with two parameters J and K as 
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where cr^ denote Pauli matrices. The last term means 
the Zeeman energy, where fi denotes the nuclear magnetic 
moment of '^He and B the magnetic field. The parameter 
J(= J3 — J2/2) is negative for most of densities, (but it 
can change the sign,) and K(= — J4/4) is always positive 
(K > 0). The first and the second summations run over 
all pairs of nearest neighbors and all minimum diamond 
clusters, respectively. The explicit form of hp for four 
sites (1,2,3,4) is 

/lp = 4(P4 + iT')-l 

= ^ (Ti ■ (Tj + {(Tl ■ Cr2){(T3 ■ 0-4) 

l<i<j<4 

-I- (cTi • (T4){cr2 ■ 0-3) - (CTI • (T3)(cr2 • (T4), (2) 

where (1,3) and (2, 4) are cj^^onal bonds of the diamond. 



The WKB approximationsQc show that exchange param- 
eters vary depending on the particle density: At high 
densities, the exchange J(< 0) is dominant, which mainly 
originates from the three spin exchange. As lowering the 
density, the ratio \K/J\ increases rapidly and hence the 
four spin exchange K{> 0) becomes important. This 
density dependence is-.caasistent with experimental re- 
sults of susceptibility.Ert£l Since multiple-spin exchanges 
produce frustration by themselves and strong copipeti- 
tions between exchanges also introduce frustrationjlj this 
model is expected to show various complex magnetic be- 
haviors. 

In a previous paper, we studied the ground state of this 
model in the classical limit and found various phases :E3 
(a) For J < —8K, the ground state shows the perfect fer- 
romagnetism. (b) For —8K < J < —8K/3, ground states 
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are highly degenerate. This degeneracy is a non-trivial 
one. (c) In — 8i^/3 < J < 25K/3, the ground state has a 
four-sublattice structure with zerojaiagnetization, which 
we call as the tetrahedral structure.^ (d) For 25K/3 < J, 
the ground state is the so-called 120° structure. Thus 
novel phases (b) and (c) appear due to the four-spin ex- 
change interactions. In the region (c), we predicted chiral 
symmetry breaking at a finite temperaturecj 

The intermediate phase (b) seems to correspond to the 
parameter region of 2D low-density solid "^He. The pa- 
rameters J and K in the low density region are estimated 
as J ~ — 1.5(mK) and ~ 0.2(mK) from susceptibility 
and specific-heat data^ and J/I^ —4 ± 2 from path 
integral Monte Carlo simulationsET t3. The five- and six- 
spin exchanges are also estimated to be comparable with 
four-spin one. This parameter region almost belongs to 
the phase (b) pf our study in the classical limit. Our 
previous studyllj shows that competitions between the 
two- and four-spin exchanges are strong in this phase 
and many kinds of ground states exist due to frustration. 
Furthermore, we found that under the magnetic field a 
plateau appears at m/mgat = 1/2 in the magnetization 
curve, where rrisat denotes the saturated magnetization. 
These various unusual phenomena occur due to frustra- 
tion caused by the multi-spin exchanges. In this paper we 
study this phase further, especially considering quantum 
effects. Finite-temperature effects are also discussed. 

In section II we summarize the results for the phase (b) 
of the classical model. Among degenerate ground states, 
one coUinear state, which we call uuud state, has the 
largest magnetization, m/mgat = 1/2, and other states 
show m/TOsat < 1/2. In section III we discuss the quan- 
tum model. Due to quantum effects, the uuud state dis- 
appears from the ground state at i? = and the ground 
state belongs to the S — Q space. Magnetization process 
of the quantum model is studied in section IV. Under 
the magnetic field, the uuud ordered state becomes sta- 
ble, since it has the largest magnetization, and makes 
a plateau at m/mgat = 1/2 in the magnetization curve. 
We discuss thermal effects in section V. Under the mag- 
netic field, which realizes the uuud ground state, the sys- 
tem shows a finite-temperature phase transition due to 
breakdown of translational symmetry. Section VI con- 
tains summaries and discussions. 



II. CLASSICAL LIMIT 

We studied|-Lhe ground state of the model (|l|) in the 
classical limit,Ej where the Pauli matrices are replaced to 
unit vectors (uf, uf, uf) with |uf p -|- jwj'p -|- jwf p = 1. 
We searched the ground state using the mean-field theory 
and studied finite-size systems with Monte Carlo method. 
We searched the ground state restricting ourselves to 
spin configurations with up to four-sublattice structures 
within the mean-field theory. To take into account larger 
sublattice structures, we studied larger finite-size sys- 




FIG. 1. Spin configuration of (a) uuud state with four sub- 
lattices and (b) a coplanar state with nine sublattices. 

tems and searched the minimum energy state with Monte 
Carlo method, gradually decreasing temperatures. Here 
we only discuss the phase (b), where the parameters are 
in SK < J < —8K/3 and competitions of the two- and 
four-spin exchanges are strong. 

If there is no magnetic field, ground states are highly 
degenerate for —8K < J < —8K/3. We list some ground 
states which we found: 

1. A collinear state with a four-sublattice structure. 
Up spins are on three sublattices and down spin 
on the other (see Fig. |l|a). We call this state uuud 
state. 

2. A coplanar state with nine sublattices, whose spin 
configuration is constructed from three kinds of 
spin vectors (see Fig. |l|b). 

3. Many other ground states can be made out of the 
uuud state by reversing all up spins to down on 
some parallel straight lines that consist of only up 
spins (see Fig. ||). 

All these states have the minimum energy E/N = —3K. 
The uuud state has the largest magnetization m/mgat = 
1/2 among these states. The coplanar state shows = 
and the third series of ground states have a variety 
of magnetization between — 1/2 < m/msat < 1/2. The 
number of degeneracy except for the SO(3) symmetry is 
at least of order 2^/^ on the LxL triangular lattice, which 
comes from line degrees of freedom for spin fiips. The 
number of the states in the = ML sector, where M is 
an integer in —L/2 < M < L/2, is at least l/2C(l-2M)/a 
and thus the degeneracy is largest in the = Q sector. 
This degeneracy does not originate from the symmetry 
of the Hamiltonian and instead comes from frustration 
effects. We suspect that strong frustration makes density 
of states large near the lowest energy and hence non- 
trivial degeneracy appears in the ground states. 

Another characteristic property also appears in excita- 
tions. For the ground states in group 3, a huge number 
of excitations have extremely low energy, which is very 
close to the ground state energy. Consider a stripe exci- 
tation which is shown on Fig. g|. On the ith line, rotate 
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FIG. 2. Spin configuration of two kinds of ground states. 
Up (down) spins are denoted by white (black) circles. The 
upper configuration is that of the uuud state. The lower one 
is obtained by reversing the up spins to down on the dashed 
lines of the upper one. 
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FIG. 3. A stripe excitation of the uuud state. Spins on the 
ith line are rotated with angle 6, and on the i — 1th and i + lth 
fines with angle 6/2. 
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FIG. 4. Parameter dependence of phase boundaries of the 
m/irinai, — 1/2 state in the magnetization process at T = 0, 
where the 1/2-plateau appears in the region surrounded by 
data. The solid line denotes the result from the mean-field 
theory in the classical limit and other data with lines denote 
the results of the quantum model {S — 1/2) on finite-size 
systems. 

spins with angle 9 and, on the i — 1th and i + 1th lines, 
rotate spins with angle 0/2, where ith line contains only 
up spins. For small 9, expanding the excitation energy 
with 9, one can find that, 0{9^) terms of the excitation 
energy vanish and the leading term starts from the order 
0{9'^), if all spins on both i — 2th and i + 2th lines di- 
rect upward. We note that line excitation energy usually 
depends on the angle 9 in the quadratic form. The condi- 
tion for the low excitations of this kind to appear is that 
all spins on three lines in next neighbors, i.e. on i — 2th, 
zth and i + 2th lines, are in the same direction with each 
other. We hence find that these low energy excitations 
exist for most of ground states in group 3 and there are 
a huge number of low-lying excitations of this kind. 

By applying a weak magnetic field, this degeneracy 
quickly disappears. Since the uuud state has the largest 
magnetization among the degenerate states, the uuud 
state is stable under the magnetic field and it has the low- 
est energy. Figure ^ shows the parameter region where 
the uuud state becomes the ground state. The uuud state 
remains to be the ground state up to finite magnitude of 
the magnetic field, which means the magnetic suscepti- 
bility to be vanishing in this phase. This behavior occurs 
due to singularity of coUinear states. Since all spin vec- 
tors in the uuud state are parallel to the magnetic field, 
spins can be rigid against the field. The uuud state hence 
makes a plateau at the half magnetization m/rrisat = 1/2 
in the magnetization process (see Fig. ||). These results 
are obtained with the mean-field theory and also con- 
firmed with Monte Carlo method. The uuud state can 
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FIG. 5. Magnetization process at T = for J = —4 and 
K — 1. The dotted line denotes the result from the mean- field 
theory in the classical limittJ and other lines denote the re- 
sults of the quantum model {S — 1/2) on finite-size systems. 



stably exist even at low but finite temperatures and it 
also shows a finite-temperature phase transition due to 
the breakdown of the translational symmetry. We will 
further discuss finite-temperature properties in section 
V. 



III. QUANTUM (S = 1/2) MODEL WITH B = 0. 

When there are various degenerate ground states in 
the classical limit, quantum effects play an essential role 
in forming the ground state of the quantum model. One 
possibility is that a new quantum ground state appears 
due to tunneling between the classical states. This mix- 
ture of states can occur between the degenerate ground 
states in the same S*^ sector. This possibility was dis- 
cussed as an origin of the disordered-state which is ob- 
served in low-density solid ^He films JIj'EI Another possi- 
bility is that one of the degenerate ground states may be 
selected due to quantum effects. 

To test these possibilities, we first examine the uuud 
state, which is one of the classical ground states, with the 
spin-wave approximation. Using the Holstein-Primakoff 
transformation, we expand the Hamiltonian up to the 
quadratic form of bosons 



n = -3NK + J2 Ai'DkAk 

k 

with aI = (al,bl,cl,d^k) and 
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FIG. 6. The lowest mode of spin-wave spectrum on the 
uuud ordered state. The contours of constant energy are writ- 
ten on the bottom plane. Zero modes appear on the three 
lines, kx = and ky 
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where a, 5, c and d denote bosons for four sublattices and 

Br,r' = 4:{J + 2K)cosk- r + SXcos fc • (r - r'), 
Cr ^ 4:{J + SK)cosk-r. (5) 

Numerioaily diagonalizing this Hamiltonian with Copla's 
method,^ we evaluate the excitation spectrum of spin 
waves.E3 The spectrum has four branches of excitation 
modes in the 1st Brillouin zone of the four-sublattice 
structure. The lowest mode shows an ill behavior (see 
Fig. H); it has flat modes along three lines in momen- 
tum space, kx = and ky = ±\/3kx/2, and furthermore 
the dispersion curves across the three lines have a form 
\k — ko\^. These modes correspond to the line excitations 
that we have shown for the classical model in Section II. 
This flat mode suggests that the uuud state does not 
have spin stiffness. Higher-order terms of spin-wave ex- 
pansions may destroy the uuud state due to non-linear 
effects. 

We also studied the ground state of finite-size systems 
with exact-diagonalization method. The systems with 
the size N = 12, 16, 20, 24 and 28 are treated. (See 
Appendix A for shapes of the finite-size clusters.) The 
results reveal that the ground state belongs to the 5 = 
space and hence it is not the uuud state for —7K < J, 
which almost covers thfi_phase (b). Recent numerical 
study by Misguich et alo also suggests that the ground 
state in the same parameter region is spin liquid with 
5 = 0. 
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FIG. 7. Size dependence of the lower- and upper-critical 
fields of the magnetization plateau at Sfotai ~ for the 

model with J = — 4 and K = 1. 



All the above results are consistent with each other 
and give a unique picture. The ground state belongs to 
the 5 = space and the uuud state has a little higher 
energy than the ground state due to quantum effects. 
The uuud state thus disappears from the ground state in 
the quantum model. But it again becomes stable under 
the magnetic field. We will discuss this point in the next 
section. 



IV. QUANTUM (S = 1/2) MODEL UNDER THE 
MAGNETIC FIELD (B > 0). 

To test the appearance of magnetization plateau at 
m/msat = 1/2 which we found in the classical limit, we 
investigate ground states of the quantum model ([|) in 
finite-size systems under the magnetic field. We study 
finite-size (N < 28) systems with periodic-boundary con- 
ditions. Figure ^ shows the magnetization process of 
finite-size systems with J = —4 and K ~ I. Though 
the magnetization increases stepwisely due to finite-size 
effects, there clearly exists a broad plateau at m/msat = 
1/2 in every-size data. Width between the lower and 
upper critical fields of the plateau does not vanish, as 
the system size increases, and remains to be significantly 
large (see Fig. |^) . This result strongly suggests that this 
magnetization plateau survives in the thermodynamic 
limit. To examine whether this m/msat = 1/2 state has 
the ituttd long-range order, we consider the following uuud 
order parameter 



^ieA ieB iec ieo 
and calculate long-range order 



(6) 



FIG. 8. Size dependence of long-range order of the uuud 
structure, {{O^) -~ {Of)/N'^, in the ground state of the 
'S'totai = -^/4 space. Black circles denote data for the model 
with J — —4 and K = 1. Data for the Heisenberg antiferro- 
magnet {J = 1, K = 0) are also shown with triangle symbols 
for a comparison. 
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in the ground state of the S^^^^y — N/4 space. Re- 
sults are shown in Fig. which clearly suggest that 
data are extrapolated to a finite value in the N ^ oo 
limit for J = —AK. The extrapolated value is estimated 
as about 0.03 in the N oo limit. Long-range order 
of the uuud structure brings breakdown of translational 
symmetry in the thermodynamic limit. A careful esti- 
mation (see Appendix B) leads that, when the trans- 
lational symmetry is spontaneous broken, the expecta- 
tion value of spin on each sublattice is (a'f/2) = 0.45 
for i G A-, B- or C-sublattice, and ((Tf/2) = —0.35 for 
i e Z3-sublattice, respectively. The ground state in the 
•S'totai = space thus has a rigid uuud long-range or- 
der and deviation of the sublattice magnetization from 
the classical value is small. By the way, the total mag- 
netization in the uuud ordered state does not change 
from the classical value m/mgat = 1/2 by quantum ef- 
fects. Note that the expectation values of spin satisfy 
^;(crf/2)/Ar = (3x0.45-G.35)/4= 1/4. The spin-wave 
analysis also does not give any quantum correction to the 
total magnetization in the uuud state. 

Magnetization plateaus of two-dimensional systems 
have been observed both theoretically and experimen- 
tally. On the triangular lattice, a magnetization plateau 
was observed at m/ni—t — 1/3 in the measurements of 
CgEucJ and CsCuCl3,E3 and in theoreti&aJ studies of anti- 
ferromagnets (with two-spin exchange) .B3~Ea The plateau 
comes from the three-sublattice uud ordered state. The 
quantum correctioDL-tp the magnetization m/mgat = 1/3 
is also vanishing. E3'E3 For MSE models, a similar mag- 
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netic plateau was also observed at m/msat = 1/2 by 
finite-size studies. Roger and Hetherington first discov- 
ered the magnetic plateau in ajMSE model with four-spin 
exchange on the square latticccjEO (In the mean-field ap- 
proximation, the ground state of this plateau is also the 
uuud state. (See Ref. ^)) Recently Misguich et al. also 
observed a plateau in a model with two-, four-, five- and 
six-spin exchanges on the triangular lattice. We believe 
that these plateaus in MSE models come from the uuud 
ordered state and that the four-spin exchange is the main 
cause for the plateau as well. 

Next, we study the excitations in the •S'^ot^i — 

N/A 

space. Figure ^ shows excitation energy of up to fourth 
excited state. Three low-lying excited states have energy 
very close to the ground-state one and they converge to 
the ground state as the system size is enlarged. Above 
them there is a large gap, which seems not to vanish in 
the N oo limit. Thus four states exist around the 
lowest level and are clearly separated from other excited 
states. The number of low-lying levels, four, is equal 
to the degeneracy of uuud states, which comes from the 
spatial translation of sublattices. Furthermore, the four 
low-lying states have the same translation group as the 
following schematic states, respectively, 

|i) = (l + 7^l + 7^2 + 7^l7^2)| ttti), 
|2) = (l + 7^l-7^2-7^l7^2)|TTTi), 
|3) = (l-7^l + 7^2-7^l7^2)| ttti), 
|4) = (l-7^l -7^2 + 7^l7^2)| TTTi). (8) 



where | ttti) denotes one of the uuud states, 
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)i SSjen I Djy and TZi {i = 1, 2) mean the translation of 
sites by unit vectors . We hence believe that these four 
levels form four ground states which have uuud order in 
the thermodynamic limit and translational invariance is 
spontaneously broken in the ground states. This argu- 
ment leads to a conclusion that the fifth low-lying state 
corresponds to the lowest excitation in the thermody- 
namic limit and hence the spin excitation spectrum has 
a finite gap in the same S^^^^y sector. 

The appearance of the magnetization plateau and the 
excitation gap can be understood by introducing a parti- 
cle picture into the spin model.Ell Let us regard the down 
spins as particles moving in background of up spins and 
the magnetic field u^B as minus of the chemical poten- 
tial of the particle.E3 Then the present spin system is 
mapped to a hard-core boson system and we can recog- 
nize the uuud ordered state as a Mott insulating state 
with charge-density wave (CDW). Note that this den- 
sity wave can be ordered by repulsion between particles 
on nearest- or next-nearest-neighbor pairs of sites, which 
originates from the four-spin exchange interaction. In in- 
compressible CDW, particles are insulating and charge- 
density excitations have a finite gap. And the compress- 
ibility is vanishing. Through the mapping, these features 
correspond to the finite excitation gap and the magneti- 
zation plateau of the original spin model. Thus the mag- 
netization plateau can be understood as an insulating- 
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FIG. 9. Size dependence of excitation energy of four 
low- lying excited states in the S'totai ~ N/A space for the 
model with J = —4 and K = 1. The symbol for the third 
excited state almost overlaps those for the first and second 
states. 



conducting transition in the particle picture. We hence 
have an explanation for the plateau from the particle 
limit (5 = 1/2). As we mentioned in section II, the 
appearance of the plateau in the 5 — > oo limit originates 
from the rigidity of the coUinear state. If the scenario in 
the particle picture is valid, the mechanism for the ap- 
pearance of plateau is understood both from the classical 
limit {S oo) and the particle limit {S ~ 1/2). Note 
that, in both cases, the uuud order is the key property 
for the appearance of the magnetization plateau. 

The 1/2-plateau appears for a wide parameter region. 
In Fig. 0, we show the phase boundary of the parame- 
ter region of the plateau. Except for weak magnetic field 
cases, the region of the m/mgat = 1/2 phase becomes 
wider in comparison with the classical model. Quantum 
effects thus stabilize the plateau and enhance its appear- 
ance. 



V. THERMAL EFFECTS 

It is also important to discuss thermal effects to the 
magnetization plateau for a comparison with experiments 
at finite temperatures. A favorable property of the uuud 
order is that it accompanies a phase transition at a finite 
temperature. Since the uuud order is a discrete symmetry 
breaking, i.e., the number of ground states is four under 
the magnetic field, the symmetry breakdown occurs even 
at low but finite temperatures in two dimensions. The 
uuud order hence survives against thermal fluctuations 
and a phase transition occurs at a finite temperature. 
(Note that this argument does not depend on whether 
the system is quantum or classical.) Fortunately, the 
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FIG. 10. Temperature dependence of the specific lieat of 
tlie classical model with J = —4, K = 1 and jj,B = 10. Monte 
Carlo simulations were done on triangular lattices with finite 
size iV = 12L^ 



magnetization plateau survives even at low temperatures 
due to ordering of the uuud structure. Magnetization in 
the uuud ordered state is close to m/msat = 1/2 at low 
temperatures and it makes a plateau (or a shoulder) near 
the half of the saturated magnetization. This argument 
comes from only the degeneracy of the ground states, 
using the aspect of universality in phase transition. 

To confirm the above argument, we study the finite- 
temperature properties in the classical limit using Monte 
Carlo simulations. We believe that critical properties 
of the phase transition at finite temperatures are gov- 
erned by thermal fluctuations and quantum effects do 
not change its universality, though the value of the criti- 
cal temperature or order parameter will deviate from the 
quantum one. Monte Carlo simulations were performed 
with the Metropolis algorithm. If a spin flip is rejected, 
we randomly rotate the spin about the local molecular 
field. We construct finite-size systems with a unit cluster 
which has 12 sites. The system sizes are N = 12L^ with 
L — A,6, 8, 12 and 16 with periodic-boundary conditions. 
After discarding initial 30000-50000 Monte Carlo steps 
per spin (MCS) for equihbration, subsequent 3 x 10^- 
5 X 10^ MCS are used to calculate the average. 

Monte Carlo simulations were done for the model un- 
der a magnetic field which realizes the uuud grovLnd state. 
Figure |l^ shows specific-heat data for the case J = —4, 
K = 1 and fiB — 10. The data show a sharp diver- 
gence around T — 1.9K. We also found that the uuud 
long-range order exists below the critical temperature. 
Near the critical temperature, the magnetization curve 
is still rounded and smooth. As lowering the tempera- 
ture, the slope of the curve decreases and becomes flat 
around, m/msat = 1/2, as shown in fig. |l^. Misguich 
et al.E3 also demonstrated stability of the plateau at fi- 




T = K 
T=K/2 
T=K/4 
T = (mean field) 



10 15 20 25 30 35 40 



I^B/K 



FIG. 11. Magnetization process of the classical model with 
J = —4 and K — 1 ai finite temperatures T/K — 0.25, 
0.5 and 1.0. The magnetization process at zepei temperature, 
which we obtained with the mean-field theoryjij is also shown 
for a comparison. 



nite low temperatures in the quantum MSB model on a 
finite-size (TV = 16) lattice. Similarly, in the Heisenberg 
antiferromagnet on the triangular lattice, the magnetiza- 
tion plateau at m/rTigat = 1/3 was successfully observed 
at finite temjjeratures in the measurements of jCeEuO 
and CsCuCIsJHj and in Monte Carlo simulationsE3 

The critical temperature is estimated as Tc = 1.9K for 
the model with J = -4K{< 0) and fiB = lOK. For 
a weak magnetic field = 5K, the estimate of Tc is 
about 1.7 K. These values can change due to quantum 
effects in the = 1/2 model. From the degeneracy of the 
ground states, the phase transition is expected to be of 
second order and to belong to the 4-state Potts universal- 
ity class. The finite-size scaling analysis, however, shows 
that critical exponent a is much larger than thje expected 
value a = 2/3 for the 4-state Potts model.a Moreover 
distribution of energy histogram in Monte Carlo simu- 
lations has two peaks at the critical temperature, which 
suggests that the phase transition is of first order. A 
similar deviation of the universality was also found in 
the phase transition of iChital symmetry breaking in the 
model ([|) for strong A'.EJ^L3 These deviations may come 
from frustration effects, or singularity of four-body in- 
teractions. Critical properties of these phase transitions 
will be discussed further in a forthcoming paper. 

The discrete symmetry of the uuud order in a mag- 
netic field accompanies two favorable properties that a 
sharp phase transition occurs at a finite temperature and 
the magnetization plateau stably exists at low but finite 
temperatures. We expect that these properties make the 
experimental observation of the magnetization plateau 
possible. 
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VI. SUMMARY AND DISCUSSION 
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In this paper, we examined the appearance of the mag- 
netization plateau at m/m-sat = 1/2 in a 2D MSE model 
on the triangular lattice, which we predicted in our previ- 
ous work. This plateau appears when the two- and four- 
spin exchange interactions compete strongly, which may 
be realized in 2D low-density solid '^He, and if the five- 
and six-exchanges are not too strong. The four-spin ex- 
change is important and relevant to make the plateau at 
m/msat = 1/2. Since ordering of Mwitd structure is break- 
down of discrete symmetry, it accompanies a phase tran- 
sition at a finite temperature. Then the plateau would 
appear below the critical temperature. Experimental ob- 
servation of this plateau in the magnetization process of 
the solid '^He films will confirm that the four-spin ex- 
change is strong and important. 

Finally we discuss possibility of observing plateau in 
experiments. Comparing our results with measurements 
of susceptibility and specific heat, we estimate the effec- 
tive exchanges J and K for low-density films. Suscep- 
tibility measurements at low densities show that the ef- 
fective coupling Jy. is antiferromagnetic, > 0, though 
it is close to vanishing .□ In our model (|^), the effective 
coupling behaves as = 2{J + 6K) and hence we can 
estimate as J ^ —6K. Specific heat measurements indi- 
cate that energy (spin) gap is small or vanishing,El where 
the system might be close to the critical point -between 
the ferromagnetic phase and the liquid phase.cj Since 
the critical point is around J/ K ~ —8 in the model (|^), 
the exchange parameters can be estimated as J ^ — 
from the specific-heat data. The value J/ K seems to be 
close to the boundary of the region where 1/2-plateau 
appears (see Fig. |4|), though two estimates are not con- 
sistent with each other. In real systems, more-than-four- 
spin exchanges may move the boundary to right or left 
(in Fig. For quantitative arguments, further studies 
will be needed on the effect of five- or six-spin exchanges. 
To make the plateau stable and wider, we would need 
lower-density solid films and enlarge the four-spin ex- 
change effect. It was reported that, by preplating HD on 
grafoil, '^He atops-solidify at lower density than double 
layer '^He films.t^Hla Thus solid ^He preplated HD may 
be a plausible candidate for observing the magnetization 
plateau. In our model with J — —AK{< 0), the lower 
critical field of the plateau is estimated about AK/ fi. Set- 
ting the parameter a,s K = l.O(mK) and /i = 2.13/iN, 
where //n — 0.366(mK/T), we estimate the lower criti- 
cal field as Be — 5[T]. (The value may be changed by 
other multiple-spin exchanges, i.e. six-spin exchange, in 
real systems.) Since this magnitude of the field is accessi- 
ble with the present experimental equipments, we expect 
experimental verification of the magnetization process to 
be possible in 2D low-density solid '^He. 
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APPENDIX A: CLUSTERS USED IN 
FINITE-SIZE STUDIES. 

In the numerical study of the S = 1/2 model on finite- 
size systems, we used finite-size clusters with the periodic 
boundary condition which matches the four-sublattice 
structure. Clusters are set on the triangular lattice as 
shown in Fig. n2. 




FIG. 12. Finite-size clusters which we used in the exact 
diagonalization study. 



APPENDIX B: EVALUATION OF 
EXPECTATION VALUES OF THE SUBLATTICE 
MAGNETIZATION. 

Estimation of the sublattice magnetization from the 
long-range order (^ needs a careful consideration about 
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symmetry breaking. In the thermodynamic hmit, trans- 
lational symmetry is spontaneously broken in a pure nat- 
ural ground state with uuud order. But the ground states 
of finite-size systems become mixed symmetric ones be- 
cause of finite-size effects. 

Here we write the mixed symmetric ground state in the 
'S'totai = space as We also consider a state in 
which translational symmetry is broken and whose ther- 
modynamic limit is the natural pure ground state. By 
the space translation of the state, \ipi) relates to other 
three uuud ordered states in the form 

1^2) = 1^3) = 7^2lV'l), IV'4) = 7^l7^2|V'l), 

(Bl) 

where TZi {Ti-2) denotes spatial translation by ei (62). In 
the mixed state, there is translational invariance. 



\Sf 



for any site i. 



(B2) 



On the other hand, in the uuud ordered state \'ipi), the 
translational symmetry is broken, 



hm 10^1^2) 



In the same way, 



lim {^,\0'\^,) 



hm -(V2|0|V'2) 

N—^oo iV 



1 



-{mi + TO2) 



hm -^{^,\0^\^p^) 

N—^oo J\J 



(B9) 



-(mi -I- 7712) 



(BIO) 



Inserting these relations into eq. (B6) and using eq. 
we have 



48mi 



24toi + 1.08 = 



(Bll) 



and then we obtain two solutions (mi, 7712) — 
(0.45,-0.35) and (0.005,0.985). From the constraints 
|mi| < 0.5 and |m2| < 0.5, we conclude that the former 
one is the physical solution. 



{ipi\S-\ipi) = mi 

(V'll'S'f IV'i) = ^2 



for ie A,B or C, (B3) 
iorieD, (B4) 



where mi ^ m2. Since the total magnetization is N/A, 
the sublattice magnetization mi and m2 satisfy 



-mi + -m2 = -. 



(B5) 



It is expected that the symmetric mixed state can be 
decomposed into the pure states in the thermodynamic 
limit in the form = (ji/'i) + e*^|?/'2) + e^'l'l^s) + 
e^'^ lip 4))/ 2, where 9, ip and (j) denote arbitrary real num- 
bers. The following relation hence holds in the large N 
limit: 



^(0|O2|0) = ^((^i|O2|^i) 



(B6) 



From the numerical calculations in section IV, we have 



hm ij(0|e)2|0)^O.O3 
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(B7) 



For the pure state, the clustering property of the state 
leads 



N 



hm -L(v,i|02|^^) = ( ^(V'lPl^i; 



N 



-(3m,i — m2) 



(B8) 
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